Fractional quantum Hall effect in graphene 
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Unlike regular electron spin, the pseudospin degeneracy of Fermi points in graphene does not 
couple directly to magnetic field. Therefore, graphene provides a natural vehicle to observe the 
integral and fractional quantum Hall physics in an elusive limit analogous to zero Zeeman split- 
ting in GaAs systems. This limit can exhibit new integral plateaus arising from interactions, 
large pseudoskyrmions, fractional sequences, even/odd numerator effects, composite-fermion pseu- 
doskyrmions, and a pseudospin-singlet composite-fermion Fermi sea. The Dirac nature of the _B = 
spectrum, which induces qualitative changes in the overall spectrum, has no bearing on the frac- 
tional quantum Hall efi^ect in the n = Landau level of graphene. The second Landau level of 
graphene is predicted to show more robust fractional quantum Hall effect than the second Landau 
level of GaAs. 

PACS numbers: 73.63.Fg, 73.23.-b, 81.05.Uw 



The discovery of integral quantum Hall plateaus in 
graphene poses several questions: Does the fractional 
quantum Hall effect (FQHE) also occur in graphene, and, 
if so, what is its character? Does it resemble the in- 
tegral quantum Hall effect (IQHE) in graphene? How 
does it compare to the FQHE in GaAs? What other 
physics does interaction produce? Even if the Zeeman 
energy is large enough to fully polarize the (real) spin 
of the low energy state, as assumed throughout this ar- 
ticle, the graphene Fermi point degeneracy, which takes 
the form of a pseudospin 1/2 character, remains. Strik- 
ingly, despite the qualitative changes in the Landau level 
(LL) spectrum that arise from the Dirac dispersion of the 
low energy states at zero field, the physics of the n = 
graphene Landau level is identical to that of electrons 
in the n — Q GaAs Landau level, but with the Zeeman 
splitting set identically to zero; the pseudospin, which 
is also a sublattice index for n = 0, plays the role of 
the traditional electron spin. This leads to the predic- 
tion of the formation of composite fermions, FQHE at 
v — n/(2pn± 1), odd versus even numerator effects, and 
a pseudospin singlet composite fermion Fermi sea. Gi- 
ant pseudoskyrmions are predicted for the ferromagnetic 
integral and fractional quantum Hall states. Our numer- 
ical diagonalizations show that, interestingly, the FQHE 
of graphene in the n — I Landau level is more robust 
than that in the n = 1 GaAs Landau level, due to more 
favorable interaction pseudopotentials. We assume below 
that the pseudospin symmetry is exact, and also that LL 
mixing is negligible. 

In graphene, the two Fermi points, each with a two- 
fold band degeneracy, can be described by a low-energy 
continuum approximation with a four-component enve- 
lope wave function whose components are labelled by 
a Fermi-point pseudospin = ±1 and a sublattice. The 
Hamiltonian for one pseudospin component is 0,01 
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where vp is the Fermi velocity and II = p + {e/c)A. 
Letting z = x — iy, choosing the symmetric gauge A — 
{—By/2,Bx/2), and introducing the Landau level rais- 
ing and lowering operators, ~ (^/2 ~ 2dz) /V^ and 
a = {z/2 + 2dz) / a/2, the Hamiltonian becomes 
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The eigenvalue problem is conveniently formulated in 
terms of , which can be represented as LL number 
operators: 
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Any 5* = ( ^'"^ I is an eigenvector of tP' , where 



I3'nn,m2 

'>ln,m are the standard LL eigenfunctions in GaAs (n = 
0, 1, • • • is the LL index, and m is the angular momen- 
tum index), i?*? = i?„'I' requires nii = TO2 and fixes 
[3/ a, with the result that normalized eigenvectors take 
the form (sgn(O) = by convention) 
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Before discussing the FQHE(|, we describe how in- 
teractions are expected to affect the integral QHE in 
graphene. For noninteracting electrons, when both the 
spin and the pseudospin degeneracies are present, the 
Hall plateaus have been predicted and seen at 
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For magnetic fields large enough to resolve the spin bands 
the Hall plateaus occur at 

RH = h/{]e'), J -...,-2, -1,0, 1,2, 4, 6, 8,... (8) 

where j = implies a band insulator. However, all inte- 
gral plateaus should, in principle, become possible for 
interacting electrons. Consider a magnetic field large 
enough to lift the spin degeneracy at a filling v = 2n + l 
with n = 1,2,3,... where no QHE occurs for noninter- 
acting electrons. For interacting electrons, the rotational 
symmetry in pseudospin space is spontaneously broken 
due to pseudospin exchange, with the pseudospin magne- 
tization picking an arbitrary direction. While this system 
supports gapless pseudospin- wave excitation, the charged 
excitations have a gap, thus producing a Hall plateau. 
FQHE in the n = Landau level of graphene: The 
unusual electronic dispersion of graphene around E = 
is reflected in both the structure of Landau levels and 
in the offset of the Hall conductance staircase. Nonethe- 
less, the actual wave functions of electrons in the n = 
manifold are identical to those in the conventional low- 
est Landau level of GaAs. Consequently, the FQHE in 
the n = Landau level of graphene at a high magnetic 
field maps onto FQHE of electrons in GaAs with zero 
Zeeman energy 5] , as corroborated and extended by the 
numerical results described below. 

The basic physics in the n = graphene LL, therefore, 
is the same as in GaAs. Each electron captures an even 
number J_2p) of quantized vortices to become a composite 
fermionp. The Berry phases generated by the vortices 
effectively cancel part of the external magnetic field, so 
that the dynamics of composite fermions are governed 
by a reduced magnetic field B* = B — 2pp(j)o, where p 
is the particle density in the n = Landau level and 
(/>o = hc/e is the flux quantum. The n = Landau level 
of electrons splits into Landau-like levels of composite 
fermions, whose filling factor v* , in terms of the electron 
filling i/, is given hy ly = v* /{2pv* ± 1). 

The IQHE of composite fermions for v* = n produces 
sequences of fractions 0: 

= n/(2pn± 1) . (9) 

The origin of gap, i.e. the energy required to promote 
a composite fermion into a higher CF-Landau level, is 
different for even and odd values of n. For even n, the 
ground state is a pseudospin singlet, with n/2 CF-Landau 
levels for each component of the pseudospin occupied. 
For odd n, the ground state is partially pseudospin polar- 
ized; no QHE would occur here if the composite fermions 
did not interact, but the residual interaction between 
composite fermions opens a gap. To the extent the resid- 
ual interaction is weak, one expects fractions with even 
numerators to be more robust than those with odd nu- 
merators. The excitation energies for the GaAs FQHE in 



the zero Zeeman energy limit apply to graphene FQHE 
within the n = level. Calculated gaps to creation of 
a far-separated charged quasiparticle/quasihole pair at 
v = 1/3 and (unpolarized) v = 2/5 are 0.07 and 0.04 
e^/e^B, respectively 0; the larger gap at 1/3 indicates 
the significance of inter-CF interactions. 

So long as electrons are confined to the n = Landau 
level, they have no memory of the Dirac nature of the 
zero-field dispersion, with some surprising consequences 
for the FQHE. The CF-cyclotron energy opens up ap- 
proximately linearly with B*, as expected for composite 
fermions with a parabolic dispersion, even though the 
cyclotron energy of electrons in graphene scales anoma- 
lously with B. 

In GaAs, the sequence of FQHE states at = n/(2pn± 
1) terminates as n — s- cx) in a composite- fermion Fermi sea 
at 1/ = l/2p, where the effective magnetic field vanishes 

lioj . For zero Zeeman energy ^variational calculations 
favor the spin singlet Fermi seajllj, so graphene should 
have a pseudospin-singlet CF Fermi sea dX v = l/2p. 
The CF Fermi sea in GaAs has been successfully mod- 
eled as an ordinary Fermi sea with parabolic dispersion, 
which allows one to deduce an effective mass for compos- 
ite fermions. The same should be true of the CF Fermi 
sea in graphene, in spite of the fact that electrons in 
graphene make a Dirac sea at zero magnetic field and 
have no effective mass. The singlet nature of the CF 
Fermi sea can be ascertained through a measurement of 
the Fermi wave vector, as was accomplished in GaAs sys- 



V 


N 


2Q 


D 


Do 


^(0) 


5(0) 


GaAs 


graphene 


overlap 


















5(1) 




5(1) 




1/3 


4 


9 


145 


5 





2 





2 





2 


0.99932 




5 


12 


1106 


10 





5/2 


2 


3/2 





5/2 


0.99998 




6 


15 


11588 


50 





3 





3 





3 


0.99875 




7 


18 


109138 


290 





7/2 





7/2 





7/2 


0.99974 


2/3 


4 


5 


29 


3 











2 








0.99968 




6 


8 


500 


10 








1 


3 








0.98887 




8 


11 


11483 


91 








2 


4 








0.97482 


2/5 


4 


7 


72 


4 




















0.99248 




6 


12 


3796 


28 








1 


3 








0.93622 




8 


17 


274842 


768 











4 








0.95578 


3/5 


5 


8 


226 


5 





3/2 





5/2 





3/2 


1.00000 




8 


13 


39131 


205 





2 





4 












TABLE I: Orbital angular momentum L^"' and 
spin/pseudospin S'"' of the ground states of finite sys- 
tems on a sphere at = i , | , | and | in the n = and n = 1 
Landau levels of graphene and GaAs. D is the dimension 
of the Hilbert space in the Lz ~ Sz = Q sector; Dq is the 
dimension of the L = sector. The last column gives the 
overlaps between the n = and n = 1 graphene ground 
states. 
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tems by various geometric means jlO|. 

FQHE in the n = 1 Landau level of graphene: 

The nature of FQHE depends on the Haldane pseudopo- 
tentials. In GaAs, the FQHE is essentially restricted 
to the lowest LL: very few fractions are seen in n = 1, 
and almost none in higher LLs. The mapping between 
GaAs and graphene does not hold in higher Landau lev- 
els, so different behaviors are expected. We first evaluate 
the Coulomb matrix elements within the nth graphene 
Landau level. Write |ni, toi; 712, m2; . . . ; tt-at, ttjat) for 
the product state ry„i,mi ^ r?„2,m2 ® ■■■ ® '?njv,m„ 
and II ni, mi; n2, TO2; . . . ; un, wat)) for (g) • • • ^ 

4((n,mi;n,TO2 || V \\ n,mz;n,mi)) = 

= (n, TOi; rt, m2|y|n, TO3; n, TO4) + (10) 

(?i — 1, TOi; n, m2\V\n — 1, 7713; n, 777,4) + 

{n, mi; n — 1, m2|y|7i, n — 1, m4) + 

(tt, — 1, mi; 77, — 1, m2\V\n — 1, ma; n — \, nii). 

By conservation of angular momentum, these matrix ele- 
ments are all proportional to Smi+rn2.m3+m4:- The prob- 
lem of electrons in the 77*'' graphene LL thus formally 
maps into that of the lowest GaAs LL (with two pseu- 
dospin copies) with an effective interaction defined by the 
pseudopotentials 
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and Kri"' is the effective pseudopotential for the 77''' LL 
in GaAs. As seen in Fig. QJa), Vr 



(l)graphene 



lies between 



Vm^ and Vm \ except for m = 1. To see what FQHE 
this interaction implies, we have numerically diagonal- 
ized finite systems in a geometry where N electrons move 
on the surface of a sphere whose center holds a mag- 
netic monopole of strength Q producing a flux of 2Qhc/e 
through the surface of the sphere. We use the pseudopo- 
tentials calculated above for the disk geometry; this gives 
the exact result for very large systems and is generally a 
reasonable approximation. 

As a preliminary step, we have investigated the fully 
pseudospin-polarized sector, where larger systems (with 
up to ten particles) can be studied, and found that the 
77 = 1 graphene LL behaves very similarly to the 77 = 
LL. This suggests that the FQHE in higher graphene LLs 
may be more robust than in GaAs. We next compare 
ground states for graphene systems including the pseu- 
dospin degree of freedom (but with the real spin frozen), 
to GaAs systems with zero Zeeman splitting. Tabic U 
shows the ground state quantum numbers, orbital angu- 
lar momentum L and spin or pseudospin S, for the 77 = 



and 77 = 1 LLs at several fractions. The 77 = results 
are identical for GaAs and graphene. The near perfect 
overlaps in the last column of the table indicate that the 
FQHE in the n ~ 1 graphene LL also strongly resembles 
that in the 77 = (lowest) graphene LL. (The analogous 
overlaps are rather low for the FQHE states in the 71 = 
and 77 = 1 LLs in GaAs [Ti.) The 1/3 state is fully polar- 
ized, whereas the 2/5 and 2/3 are (pseudo)spin singlet. 
At 3/5 the spin of the ground state differs from that of 
the lowest LL for N = 8; the existence or the nature of 
FQHE at this fraction remains unclear at the moment. 

Pseudoskyrmions: In GaAs quantum wells, the excita- 
tions of the V ~ 1 state for exactly zero Zeeman splitting 
are not simple particle-hole excitations but spin textures 
called skyrmions 13, 14, 15], in which half of the spins 
are reversed. However, the skyrmion size rapidly de- 
creases with increasing Zeeman energy; experimentally, 
skyrmions typically have 3 to 5 flipped spins 16, 17J. No 
skyrmions occur at 7/ = 3, 5, . . .. CF skyrmions are be- 
lieved to be relevant near = 1/3 at very small Zeeman 
energies H 00. 

In contrast, for the state in which one of the two degen- 
erate levels of the 77 = graphene manifold is fully occu- 
pied (which produces zero Hall conductance), the excita- 
tions ought to be large pseudoskyrmions. In exact diag- 
onalization studies, we find pseudoskyrmions also in the 
77 = 1, 2 graphene LLs, where the addition of one particle 
or hole to the fully pseudospin polarized state produces 
a pseudospin singlet state. No such behavior is seen for 
77 > 3; here the excitation is fully pseudospin-polarized 
{S — N/2) on the quasihole side and has a single pseu- 
dospin reversed {S — N/2 — 1) on the quasiparticle side. 
Fig. db) depicts the N dependence of the gap to creat- 
ing a pair of pseudoskyrmion and anti-pseudoskyrmion, 
computed by exact diagonalization. (We follow the con- 
vention of Ref. jl4| to define the pseudoskyrmion gap in 
terms of "neutral" quasiparticle/quasihole energies |21j.) 
Extrapolation to the thermodynamic limit yields A^*^^ = 
0.606(15), A^^' = 0.126(7), and A^ ^ = 0.18(1). The gap 
in the 77 = LL is consistent with -y/Tr/S (e^/e^s), half 
the energy required to create an ordinary particle-hole 
pair excitation 0| . We note that the pseudospin tex- 
ture in the 71 = graphene LL can be imaged directly by 
scanning tunneling microscopy, since an electron's pseu- 
dospin determines on which sublattice it resides. 

Fig.njc-d) shows the gaps as a function of the number 
of particles for several fractional filling factors in the 77th 
graphene LL. The lowest-energy charged excitations at 
v'^^^ = 1/3 are CF pseudoskyrmions (pseudo-spin singlet 
state), with an excitation energy of 0.017, to be compared 
to 0.096 for a (pseudospin reversed) particle-hole pair of 
composite fermions. (The latter gap is greater than in the 
77 = LL because the Vi pseudopotential is greater in the 
77 = 1 than in the 77 = LL.) The gaps at filling factors 
2/5 and 2/3 involve pseudospin reversal for composite 
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FIG. 1: (a) Pseudopotentials for n = and n = 1 Landau levels in graphene and GaAs. (b)-(d) Activation gaps at 
_ 1^1/3 and 2/5 in several graphene LLs. All energies are given in units of e^/eln, where e is the dielectric constant 
of the host semiconductor and Ib is the magnetic length. The gaps in (b) and (c) refer to the energy required to create a 
pseudoskyrmion-antiskyrmion pair. Gaps to pseudospin conserving (cons.) and pseudospin reversed (rev.) excitations are given 
in (d). 



fermions. For 2/5 the gaps, A^"^^ = 0.051(1) and A^^g - 
0.062(1), were obtained from the trial wave functions of 
the CF theory, evaluated by a procedure described in Ref. 

We summarize our principal conclusions. The lowest 
LL FQHE of graphene in the large Zeeman energy limit 
is equivalent to the lowest LL FQHE in GaAs in the 
zero Zeeman energy limit, terminating into a pseudospin- 
singlet Fermi sea at half filling. The effective interaction 
in Eq. (|1HI is shown to be more favorable to CF forma- 
tion in the |n| = 1 LL of graphene than in the n = 1 LL 
of GaAs. The gaps at = 1/3, 2/5 are calculated and 
FQHE is predicted at = 2/3 due to reverse flux at- 
tachment. In contrast to GaAs, skyrmions are predicted 
to occur at z^*^"^ = 1, 3 in the n = 0, 1, 2 LLs. 

While completing the manuscript, we learned from 
Philip Kim of the observation of new integral plateaus by 
Zhang et al. |23l| . Several recent theory preprints have ex- 
amined similar topics, including a field theory approach 
that obtains results similar to ours jl^l and derivations 
of the Coulomb pseudopotentials for graphene [2^ . Sup- 
port by the National Science Foundation under grants 
no. DMR-0240458, DMR-0305035 and ECS-0609243 is 
gratefully acknowledged. 
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